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This paper uses a palm-sized quadrotor to demonstrate the gas source 
localization (GSL) in obstacle regions. The chemical substance performs a 
dilute dispersion in the air, and it is challenging to locate the source position. 
The difficulty level of the GSL increases if the location is in the obstacle 
region. To estimate the gas source position, it indeed uses high 
computational cost systems. However, it leads to a new problem: a palm- 
sized quadrotor has a limited size and payload capacity. In this paper, a 
search based on the silkworm moth is employed. The algorithm works 
straightforwardly as the detection gas triggers its motion. A repulsion 
function is integrated into a searching method to compensate the algorithm 
for working in single surface obstacle regions. The validation data and the 
palm-sized quadrotor design are provided to convince the search 
effectiveness performance. 


This is an open access article under the CC BY-SA license. 


Corresponding Author: 


Muhamad Rausyan Fikri 

Department of Information Systems, Sampoerna University 
Jakarta, Indonesia 

Email: muhamad.fikri@sampoernauniversity.ac.id 


1. INTRODUCTION 

This paper is an extended work on gas source localization (GSL) [1]. In a study [1], the validation of 
GSL performance was conducted in an open-space environment. The robot's ability to discover the gas 
source is an attractive topic with various proposed methods [2]-[5]. Implementing the GSL is a challenging 
engineering problem due to the frequently diluted gas distribution. On the other hand, the performance of the 
GSL is helpful for some applications such as finding victims in a disaster area and locating the source of gas 
leakage in mining companies [6]-[9], forest fires [10], [11] and investigating fruit maturity [12], [13]. A 
comprehensive study of the GSL is available in [14]. The implementation of GSL shows some promising 
features; however, estimating the source position from the gas characteristics is challenging thanks to its 
complexity and turbulent distribution [15]. When it contacts the air [16], the challenge may increase 
significantly with a particular assumption that this task is performed in a disaster area with obstacles 
surrounding the system's range of interest (Rol). 

In the previous study, some methods for GSL from animal behavior, known as the term reactive 
searching, have been proposed. This technique is inspired by animal behavior, such as drosophila [17], the 
navigation strategies of mice [18], bio-inspired Levy Taxis [19], bat algorithm [20], and mating partner 
seeking behavior by the moth [1], [2], [21]-[23]. To name of few, 3D investigation is also conducted by [2], 
[5], [19], [24], [25]. However, the research focuses on the open space area, and few use quadrotors. On the 
other hand, some methods have already been proposed from non-animal behavior or active searching. In the 
study [3], a particle filter algorithm is employed to locate the gas in a quadrotor. Hutchinson ef al. [26] 
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employ an information-theoretic method to perform source localization, and [5] employ a fuzzy-based filter. 
The implementation requires high computational resources and impacts the palm-sized quadrotor's payload. 
The integration of GSL and vision systems using a camera is also reported in some works [27], [28], where it 
is said to implement the search using vision systems increases the payload significantly and reduces its 
agility. Vision implementation is limited to specific applications such as mapping and perceiving the source 
[29]-[32]. The drawbacks are the time duration needed to inspect the environment entirely for mapping and 
the strong object assumption for the vision-based source localization. On the other hand, the GSL in an 
unknown environment often encounters obstacles. In this situation, the palm-sized quadrotor with the 
efficient payload could be the option because it can maneuver freely, is agile, and has less computational 
cost, especially when the system integrates with the obstacle avoidance techniques. Then, developing the 
technique for searching while avoiding obstacle is essential. According to the previous works [1] and [33], 
using silkmoth behavior shows a promising result in searching performance. In this work, Russell et al. [34] 
and Kowadlo and Russell [35] also noted that the silkmoth exhibits a direct action when receiving the female 
moth's pheromone. 

Knowing the potential of searching and obstacle avoidance techniques can be implemented in this 
situation. Thus, this study proposes two contributions. The first is developing more efficient sensing systems 
than the previous [1]. Some modification is introduced to support the palm-sized quadrotor performing GSL 
in a confined area. The second contribution is developing a collaborative algorithm between the GSL and the 
obstacle avoidance algorithm. Here, the sensor position is a priority to conduct the investigation. An 
experimental validation using particle image velocimetry (PIV) is introduced. The avoidance algorithm is 
based on the repulsion force artificially sensed by the distance sensors. Therefore, the mathematical method 
of this idea is elucidated in this paper. 

The content of this paper is organized as follows: In section 2 focuses on the problem statements to 
have a clear goal of this research direction. The explanation of the silkmoth behavior mechanism of the GSL, 
data processing technique, and obstacle avoidance will be elucidated in section 3. In section 4, the evaluation 
of sensor position is performed using PIV. The description of the experimental design to evaluate the 
quadrotor performance in an obstacle region will be delivered in section 5. Finally, the discussion and 
conclusion of our study will be given in sections 6 and 7. 


2. THE PROPOSED METHOD 

In this section, three components fof the palm-sized quadrotor are discussed. The first is the 
command execution system (CES) design for GSL and obstacle avoidance, the second is GSL system 
workflow, and the third is sensors selection and layout. The dynamics and searching mechanism and the 
experimental design for searching validation are also discussed. The problem in this paper is formally 
introduced in the problem statement. 


2.1. Problem statements 

A small take-off weight capacity causes a significant payload constraint in a palm-sized quadrotor. 
Consequently, the energy resource is minimal, and the sensing system size and weight are minimal. In this 
case, the flight duration is short. Additionally, a palm-sized quadrotor must carry the other systems, such as a 
distance sensor for obstacle avoidance. Here, two main issues are considered to search in the obstacle region: 
A palm-sized quadrotor component design for gas localization and avoiding obstacles. In this case, the proper 
selection for the system can be decided carefully. Second, determining the method to maintain the searching 
in the obstacle region under the time constraint. Here, the modified bioinspired algorithm is introduced. 

The palm-sized quadrotor uses parrot mini-drone mambo FPV (Parrot Inc), where flight duration 
information in a palm-sized quadrotor without payload is 10 minutes. The flight duration becomes significantly 
shorter as the payload increases, around 5 to 8 minutes. Then, determining the quadrotor dimension was 
discussed in [1], [23], [36], where the dimension of a palm-sized quadrotor was approximately 130 mm. The 
reference for payload limitation is based on [1], with the maximum payload is must less than 40 g. 

The next step is to embed the modified bioinspired algorithm into the system. The modified 
bioinspired algorithm consists of searching based on silkmoth behavior and potential function. The 
algorithm's effectiveness will be validated through experimental trials. First, however, a brief description is 
included to explain how the palm-sized quadrotor can maneuver in the obstacle region. The environment for 
the searching method validation uses two designs. The first design is a blocked space where the obstacle 
placed in the front of the gas source aligns with the quadrotor position. The second design is a half-open 
space where the obstacle is placed in front of the gas source, but there is a gap for the gas flows directly 
toward the palm-sized quadrotor. Therefore, developing a GSL robot using a palm-sized quadrotor should 
satisfy the aim, which is the payload limit, and time limit using two design cases. The gas sensor position 
validation is done using PIV, and in this study, the searching focus is limited to the 2D area. 
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2.2. Command execution system and power 

The role of the CES is critical to determining the behavior of a palm-sized quadrotor. The CES 
consists of the Arduino block from Sparkfun and Intel Edison as the control board. Each CES component 
communicates and transmits the data from the sensors and decides the behavior for palm-sized quadrotor 
movements. For detail, please see Figure 1. The CES collects the data from the sensor by using Arduino and 
transfers it to the server through Intel Edison as the core part handles the data transfer via Wi-Fi 
communication. All processes to execute the command behavior are conducted automatically without human 
intervention, except the initial palm-sized quadrotor take-off, landing, turn-on, and turn-of-the-system power. 
Inside the CES, communication between Arduino and Intel Edison to pass the data uses serial and I2C 
communication depending on the type of sensor attached to the CES. 
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Figure 1. Diagram block of command center 


Another important part to support the CES is the power system battery. The battery to support both 
CES and palm-sized quadrotors is selected carefully. The quadrotor's battery uses the Morpillot Li-Po battery 
(Amazon.jp) with a 3.7 V 600 mAh and a weight of 9.07 g. CES's power is used Li-Po battery 3.7 V 380 
mAh, with a weight approximation of 10.3 g. The battery's power systems are separated to minimize the 
complexity of merging the different voltage consumption between a palm-sized quadrotor and CES. 
Typically, the flight drains more power instead of supporting the hardware system. 


2.3. Sensors selection and configuration 

In this study, the primary sensors are gas and distance sensors. Since the palm-sized quadrotor has a 
size limitation, however, a small size sensor is selected. Here, the gas sensor uses MEMs MiCS-5524 Gas 
Sensor Breakout (Adafruit), which has small dimensions of 20.0x12.7x3.1 mm with a weight of 0.9 g, the 
MiCS-5524 sensor has a performance almost equal to the natural silkmoth antennae [37]. 

To support the palm-sized quadrotor avoids any obstacle, a MEMS size VL53L0X is chosen where 
it has a dimension 21.0x18.0x2.8 mm with 1.3 g of the weight. The sensor also can detect the object from 
two meters away. Two meters range of detection can give time to a palm-size quadrotor to decide if the 
computational speed is low. Hence, both gas and distance sensors satisfy the criterion of small size and 
lightweight. Then, the gas sensor position needs validation, where the validation performs by conducting a 
PIV to study the wind flow's chemical distribution [38]. The analysis procedure is that the experiment should 
be conducted in a room without light or completely dark. The analysis uses a laser, and this condition aims to 
obtain a satisfactory result to clearly understand the gas distribution behavior, as illustrated in Figure 2. 

The PIV study suggests that the sensor position is located in the front of a palm-sized quadrotor. 
Thus, the gas sensor’s position is ideally mounted on the front of the quadrotor body, assuming the gas 
distribution approach from the front side. However, including Intel Edison (5.2 g) and Arduino block from 
Sparkfun (7.5 g), the total load is 38 gr which means, this result achieves less than the allowed weight of 40 gr. 
The list of components is shown in Table 1. 
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Velocity(UV) 
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Figure 2. Condition of gas distribution in same flow with the condition before and after affected by the 
rotation of palm-sized quadrotor propellers 


Table 1. Component detail 


Name 
Processing Unit Intel Edison 
Arduino block from Sparkfun 
Gas sensors MiCS-5524 
Distance Sensor VLS3L0X 
Palm-sized quadrotor parrot mini-drone mambo FPV (Parrot Inc), 


2.4. Palm-sized quadrotor searching mechanism 
2.4.1. Palm-sized quadrotor dynamics 

Consider a palm-sized quadrotor flight in two dimensions space x-y plane with z-axis is fixed. The 
dynamic of the quadrotor in two dimensions is given by (1): 


å = flaw a) 


where q=[x y y]? and u= [u, Uy Uy]? are the global position and local input control, 
respectively. As the q holds the global terms for the position, then the dynamics can be described in (2): 


x cosy -siny O][Ux 
Y\=|sinp cosp o| bs (2) 
wp 0 0 1] Luy 


Thus, the quadrotor model in (2) defines by (3): 


X = ux Cos — uy sin yp, 
ý =u, sin Y + uy cosy, (3) 


Y = uy. 


2.4.2. Sensor data pre-processing and searching mechanism 

The searching mechanism starts by obtaining data from the gas sensor and determining the action 
which satisfies the threshold after several data pre-processing. Now, assuming the sensor detection is 
sequential, then, s € [s,, s+] is the data from the gas sensor for the left and right sides, respectively. The gas 
detection has the potential that the noise disturbs the data reading. In order to reduce the noise, average 
means filter (AM) is employed, which define by (4): 


5(k) =E s(k — i) (4) 


where, S(k) is the output data sensor after filtering, and s(k — i) is the raw data sensor before noise 
reduction. However, the after-filtering data may give a true false detection. Then, the inverse model based on 
autoregressive exogenous (ARX) to estimate true input is introduced to avoid the false detection. This model 
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will simplify the determination of behavior patterns for the palm-sized quadrotor. The implementation of the 
ARX model to estimate the ethanol gas detection to both the left and the right gas sensors defines by (5): 


§(k) = bos(k) + bys(k — 1) — a(k — 1) — a2ŝ(k — 2). (5) 


The parameter of a4, az, bọ, baare the coefficient of estimation determines by linear least squares. 
The in (5) is applied to both left and right gas sensors. Thus, $ € [ŝ;,, +] is the estimation calculation for the 
gas sensor on the left and right sides. The calculation in (5) is compared with the gas sensor threshold 6 to 
obtain the logical result. We denote the logical result m € [0,1] by checking the value of ŝ define by (6): 


= a S(k) < ô 


0, 8k) > 8 (6) 


Then, the comparison between two sensors is needed to determine the quadrotor turning direction, which is 
set for each $, and §,.. The turning direction is elucidated in algorithm 1. The time t is set by 5 seconds. 


Algorithm 1. Algorithm of silkmoth pseudocode 
1 Input Gas Input stimuli, s + 0 

2 Output Behavioral movement 

3 Do calculation in equation (2) and (3) 

4: If ŝ& < 6 or $, < 6 then 

5: Perform surge 

6: If §, < &, then 

7 Do zigzag motion, turning to left 
8: Do time counting, t++ 

9: Ift == 5 then 


10: Do a loop motion, turning direction to the left side 
11: endif 

12: endif 

13: If § < $, then 

14: Do zigzag motion, turning to left 

15: Do time counting, t++ 

16: Ift ==5 then 

17: Do a loop motion, turning direction to the right side 
18: endif 

19: endif 

20: Endif 


2.4.3. Palm-sized quadrotor surge motion 
The silkmoth mechanism starts when §(k) < 5. The behavior of the silkworm is embedded into the 
sequential motion shown in Figure 3. Figure 3 shows the behaviors of silkmoth correspond to Algorithm 1. 


Loop Zigzag Surge 


F z -A 
loop | Zigzag | Surge 


Figure 3. Silkmoth mechanism 


The surge event occurs when the gas detection satisfies the (6). Thus, the palm-sized quadrotor will 
perform a surge with the surge speed is defined by (7); 


surge = K,m (7) 
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where, K, is the desired gain for a surge. This surge parameter will influence the speed of the quadrotor when 
it detects the gas. When the sensors are not detecting any odor cues. Then, the palm-sized quadrotor will 
continue to the following behavioral state to zigzag and looping. 


2.4.4. Palm-sized quadrotor zigzag and loop motion 

Determining the turning direction for zigzag and looping can be extracted from the information of 
the first sensor that detects the gas. For example, when the left sensor detects the gas, the turning direction 
after surge will change to zigzag, turning directly to the left side. As for the right-side sensor to detect the 
information, the turning direction will turn to the right side. A random turning direction if each sensor detects 
the gas substance simultaneously. The turning duration for zigzag is 5 seconds. Suppose the sensor lost the 
gas information while in the middle of the zigzag state, the behavioral motion will change to looping when 
the times already passed from 5 seconds. The t = 5 corresponds to the algorithm in Algorithm 1. 


2.4.5. Obstacle avoidance 

Obstacle avoidance is designed to ensure the safety of the palm-sized quadrotor when searching in 
an environment with an obstacle. The three distance sensors are attached following this configuration (front, 
left, and right) side of the quadrotor. To activate the quadrotor behaviors when the distance sensors detect 
any obstacle in their position. Let set the threshold 7 as the reference for the minimum distance to prevent the 
collision. The behavior of the quadrotor to prevent the collision is then described by (8); 


= (8) 


Fe d(k) <n 
0, dik) =n 


where, p is the logical operator, where p € [0,1]. The parameter of d(k) is the distance sensor parameter. 
Then, the speed to avoid collision with the object is defined by (9). 


avoid = —k,p (9) 


In (9), Ko is the desired gain of avoidance. Therefore, in (9) can be applied in the x-y plane. The negative 
sign means the effort of the quadrotor to escape from the collision. There will be two turning states when the 
palm-sized quadrotor starts searching based on the flowchart in Figure 4. The first is when it detects the gas 
information, and the second is when it detects the obstacle. When the palm-sized quadrotor approaches the 
obstacle, it will switch its behavior to obstacle avoidance. The switching behavior changes depending on 
which information the systems receive. 


CPT 


Obstacle | Surge 
Avoidance | 


Which 
sensor < 


threshold? , i Turn Turn Turn 
right i! Left Randomly Right 
Turn to Turn Turn to l i I I 
right randomly left l Zigzag 
$ l $ 
Wait for stimuli i! Loop 
— = = e e e e e i m m | E E eg 


Figure 4. Flowchart of GSL and obstacle avoidance in palm-sized quadrotor 


To add the obstacle avoidance technique, in the x-y axis should be stated where the proposed a 
repulsion-force based to improve the quadrotor capability which defines by (10). 
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Uy = AVOId net, 
= —KoPnet (10) 
= —Kop, — Kopy. 


where p; and p, represent the sensor on the palm-sized quadrotor's left and right-side logical operator, which 
has the same function as (9). Next, the input control in the y-axis is affected by one distance sensor and gas 
sensor to decide the GSL motion as surge. Thus, the u, will be (11), 


Uy = surge + avoidnet, 
= Ksm — KoPnet, (11) 
= Ksm — Kopp 


where pp is representing the logical operator of the distance sensor on the front side. For uy = f (t) where 
f(t) is piecewise continuous function. Therefore, the configuration of the quadrotor model will be (12), 


cosy -siny 0]|—KoPi — KoPr 
1} fin yp pi o| | Ksm — Kops (12) 
1 f(t) 


2.5. Experimental design 

In this study, the performance of GSL is validated in an obstacle environment. Therefore, two 
scenarios will be offered to validate the search performance: the first scenario is the obstacle in front of the 
gas source. The second scenario is that the obstacle is in the front of the source, but there is space to let the 
gas flow. The gas source's experimental design | and 2 is shown in Figures 5 (a) and 5(b). 


Palm-sized 
Quadrotor 


es 


(b) 


Figure 5. Experimental design: (a) obstacle blocks the source and (b) obstacle blocks the source with a 
certain gap 
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The experimental design is set to validate the GSL and the obstacle avoidance performance. In this 
experiment, the gas flows continuously and directly toward the quadrotor's location, the flow rate is 1L per 
min, the wind speed is 0.65 ms measured by an anemometer, and the palm-sized quadrotor’s height level is 
lm. The room temperature is set to 24 °C and the humidity level is 40%. After each trial, the exhaust fan is 
turned on to keep the environment less from ethanol pollutants, disturbing the sensor reading. Moreover, the 
quadrotor speed is set at 5 cm/s. The coefficient of ARX is set -0.981, 0.01653, 0.2833, -0.2706 for a4, 
Az, bo, b4, respectively. The threshold for gas detection is set at 0.033 for each left sensor and the right 
sensor. The duration of searching maximum is 80 seconds considering the palm-sized quadrotor time 
duration because of carrying the payload. 


3. RESULT AND DISCUSSION 

In the experiment, the quadrotor conducts the GSL in two environmental designs. Before elucidating 
the GSL performance, the total payload for GSL system design is achieved at a weight of 38 gr lower than 
the allowed payload weight of 40 gr, where the limit of 40 gr is the threshold for maximum take-off weight 
excludes the frame and the internal battery. The parameters elucidated in the experimental design section, the 
quadrotor runs using node.js programming language. The programming language is chosen because node.js 
provides a library to control the parrot mini-drone and the sensor straightforwardly. The first trial is 
conducting the GSL in the environmental design shown in Figure 5(a) to test the system's effectiveness. The 
result and trajectory obtained in this experiment based on Figure 5(a) configuration are shown in Figure 6, 
where searching performance is shown in Table 2. 

It is interesting to know that searching in blocked space, the success rate is 0%. Further 
investigation, a significant problem is faced in our platform: the compatibility of a searching algorithm that 
must compensate for the time limit during the mission. The palm-sized quadrotor's incompatibility can be 
caused by the gas source concentration dispersed to the position of the palm-sized quadrotor itself is weak. 
This result is thanks to the obstacle that blocks the source. This result confirms the experiment conducted in 
[39], which shows that the visualization of the gas flow reduces significantly. Another research that confirms 
the challenge of searching in the environment in Figure 6 was mentioned in [40]. In [40] result, it is known in 
the backside of the obstacle that the gas concentration is low thanks to the near-wake region. Therefore, the 
palm-sized quadrotor only repeated the surge-zigzag-loop pattern in this experiment until it depleted its battery. 


Experimental Design 1 


-500 0 500 1000 1500 2000 2500 
Distance (mm) 


Figure 6. Result in experimental design follows Figure 5(a) 


Table 2. GSL Performance with experimental design in Figure 5(a) 
Success Unsuccess 


Trials 0 10 
Average searching time Os >80 s 
Minimum searching time Os >80 s 
Maximum searching time Os >80s 
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On the other hand, the wind contributes to creating a wake around the obstacle. The wake 
phenomenon that occurs around the obstacle caused the gas substance to flock near the obstacle. As the gas 
substance flock in the rear part of the obstacle, the concentration increases not in the part of the obstacle 
headed to the palm-sized quadrotor. The investigation of gas searching continues to experimental design 2. 
The success rate in searching the gas source is 80% see Table 3 and Figure 7. Here, the success rate is 
significantly different from the experimental design 1. This success searching result (design in 
Figure 5(b) if we compare with [1] is increased wherein the previous research [1] approximately: 70% of 
success rate this thanks to the flow of the gas which its direction affected by the position of the obstacle. The 
searching success rate is also contributed by a gas distribution which is not diluted because of blocked by an 
obstacle. This result also confirms [40] result. To validate the GSL, we conduct experimentation in two 
experimental designs. This result confirms that the obstacle position design directly affects the unimpeded 
gas substance flow to the palm-sized quadrotor. The exciting part of experimental design 2 in Figure 5(b) is 
that the unsuccessful attempt occurred because the systems hardly recognize the gas information, which 
causes the searching time to exceed the limit. The algorithm of potential function to avoid the collision is 
successfully managed by the palm-sized quadrotor to avoid the crash. 


Table 3. GSL performance with experimental design in Figure 5(b) 
Success Unsuccess 


Trials 8 2 
Average Searching time 31.4s 81.25 s 
Minimum searching time 15s 30s 
Maximum searching time 70s 120 s 


Experimental Design 2 


n = 10, success = 8, failed = 2 
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N 
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-1000 
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Figure 7. GSL performance with design in Figure 5(b) 


4. CONCLUSION 

This paper demonstrates a GSL using a palm-sized quadrotor in the obstacle environment and 
verified through experimental validation. This paper has developed a 38-gr system of the palm-sized 
quadrotor, consisting of obstacle avoidance and GSL systems and below the allowed payload limitation. The 
performance of searching in the obstacle region varies in two designs where it achieves 0% and 80% with the 
explanation in the previous section. Based on the results, the algorithm of the insect behavior works 
appropriately to respond to the gas information. This result confirms that the GSL algorithm can work 
simultaneously with obstacle avoidance. This result also proves the GSL integration with obstacle avoidance 
(modified bio-inspired). According to the palm-sized quadrotor conditions, improving gas detection design 
and algorithm is required. Moreover, the limited payload size also needs consideration because the quadrotor 
design is relatively small. 

This study considers the drawbacks that occurred in Figure 6. The gas distribution is significantly 
reduced in configuration Figure 5(a). The bio-inspired algorithm employed in this study works if there is a 
gas detection. However, when the gas information is lost due to a particular condition, such as an obstacle. In 
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that case, the recommendation for future reference is that the system must be improved by adding an active 
searching algorithm or using a palm-sized quadrotor with a built-in visual system to improve navigation. 
Searching in 3D also can be considered to validate the performance of the proposed algorithm. 
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